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ABSTRACT 

A  Model  to  Estimate  the  Rux  of  Background  Particles 
Expected  to  Accumulate  on  the  Collector  Foils 
of  the  Interstellar  Gas  Experiment  on  the 
Long  Duration  Exposure  Facility  Spacecraft 

by 

Paul  A.  Gehred,  Master  of  Physics 
Utah  State  University,  1992 

Major  Professor.  Dr.  Don  L.  Lind 
Department:  Physics 

The  Interstellar  Gas  Experiment  (IGE)  on  board  the  Long  Duration  Exposure  Facility 
(LDEF)  exposed  metallic  foils  to  the  intersteUar  wind  to  collect  noble  gas  particles  from  outside  of 
the  solar  system.  The  terrestrial  helium  present  in  the  atmosphere  at  low  Earth  orbit  altitudes  was 
also  collected  on  the  foils,  along  with  the  interstellar  helium  particles. 

In  order  to  discriminate  between  the  desired  interstellar  helium  and  the  background  terrestrial 
helium  particles,  an  analysis  of  the  expected  distribution  of  atmospheric  helium  on  the  foils  was 
conducted.  This  analysis  used  the  principles  of  the  kinetic  theory  of  gases  to  simulate  the  state  of  the 
thermosphere  for  the  entire  69-month  flight  of  LDEF.  The  key  parameters  needed  to  apply  the 
kinetic  theoiy  of  gases  were  the  temperature  and  density  of  the  helium  at  LDEF  altitude.  This 
information  was  calculated  using  the  Mass  Spectrometer  and  Incoherent  Scatter  (MSIS)  model. 

The  entire  duration  of  the  flight  was  split  into  70  month-long  intervals.  For  each  of  these 
intervals,  mean  values  for  the  kinetic  temperature  and  helium  density  were  calculated.  These  values 
were  calculated  twice  per  month,  once  for  the  AM  period  and  once  for  the  PM  period,  in  local  solar 
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time.  The  average  altitude  and  orbital  speed  were  also  calculated  for  each  month.  The  integration 
program  was  fed  140  sets  of  values  of  these  parameters. 

The  density  gradient  of  the  accumulated  particles  was  calculated  using  the  Maxwellian 
velocity  distribution  together  with  the  limiting  geometry  of  the  foil  collector  housing.  The  collector 
housing  walls  defined  the  field  of  view  of  each  foil.  The  field  of  view  changed  slightly  with  the 
position  on  the  foil.  For  each  month,  the  integration  calculated  the  flux  at  five  different  positions 
across  the  center  of  each  foil.  The  limits  of  the  field  of  view  were  used  as  the  limits  of  integration 
in  three-dimensional  velocity  space. 

Another  variable  in  the  integration  was  the  probability  that  a  particle  which  struck  the  foil 
would  actually  adhere  to  it.  This  probability  was  a  function  of  the  velocity  of  the  particle.  An 
experiment  conducted  in  1983  provided  a  set  of  18  data  points.  These  data  points  were  fitted  to  a 
curve  called  the  eta  function,  which  provided  continuous  probability  values  for  any  chosen  velocity. 
Some  ambiguity  in  the  data  for  this  function  suggested  solving  the  integration  with  an 
underestimation  and  with  an  overestimation  of  the  eta  function  in  order  to  give  a  range  of 
uncertainty  for  the  possible  solutions. 

Results  of  the  calculations  depend  on  the  orientation  of  the  collector  with  respect  to  the 
orbital  velocity.  For  the  midpoint  of  foils  which  were  tilted  perpendicular  to  the  orbital  motion,  the 
expected  helium  accumulation  was  5(10)®  helium  atoms  per  square  meter  of  foil.  For  collector 
three,  which  was  dipped  24"  into  the  direction  of  orbital  motion,  the  expected  accumulation  was 
5(10)*  atoms  per  square  meter.  Collector  two  was  dipped  24"  out  of  the  direction  of  orbital  motion, 
and  was  expected  to  accumulate  only  6(10)^  atoms  per  square  meter.  Due  to  the  uncertainty  in  the 
eta  function  these  values  are  plus  or  minus  at  least  an  order  of  magnitude. 

The  mass  spectrometer  analysis  conducted  so  far  has  found  traces  of  ^He,  ^He,  ^®Ne,  and 
^Ne.  Since  MSIS  ‘jnly  provided  data  for  the  ^He,  estimated  values  of  the  other  isotopes  were  based 
on  the  relative  abundances  proportional  to  ^He.  The  result  of  this  research  will  hopefully  supply 
the  experimenters  with  values  to  be  subtracted  out  of  the  final  concentration.  The  remainder  will 
provide  information  on  the  density  of  the  interstellar  wind  and  the  isotope  ratios  of  the  noble  gases 


xi 

interstellar  wind.  This  information  will  aid  research  of  nucleosynthesis  and  provide  new  clues 
regarding  the  origin  of  the  universe. 

(88  pages) 


CHAPTER  1 


INTRODUCTION 

1.1  The  Interstellar  Gas  Experiment 

The  purpose  of  the  Interstellar  Gas  Experiment  (IGE)  was  to  collect  and  isotopicaliy  analyze 
the  noble  gas  compon'^nt  of  the  interstellar  medium.  The  scientific  significance  of  the  data  to  be 
derived  from  the  direct  detection  of  interstellar  particles  is  in  understanding  nucleosynthesis.  Current 
techniques  for  determining  relative  abundances  of  stellar  elements  rely  on  stellar  spectroscopy.  This 
technique  falters  when  applied  to  noble  gases  because  the  interpretation  depends  strongly  on 
assumptions  about  the  atmospheric  structure  of  the  source  [Lind  et  oL,  1988].  The  isotopic  ratios  can 
also  be  determined  more  readily  with  the  mass  spectrometer  analysis. 

Tlie  collection  device  used  in  the  JOB  was  made  up  of  seven  metallic  foils.  High  energy  noble 
gas  particles  impacting  the  foils  became  lodged  in  the  metal  lattice.  These  seven  collectors  were 
placed  on  the  Long  Duration  Exposure  Facility  (LDEF)  with  different  fields  of  view.  These  view 
angles  optimized  the  chance  for  the  authors  of  the  ^qreriment  to  recognize  the  angular  distribution 
of  the  interstellar  gas  particles.  The  solar  gravitation  field  focuses  the  interstellar  wind  into  a  type 
of  stream.  As  the  LDEF  spacecraft  revolved  into  the  leeside  of  the  sun,  it  bathed  in  this  stream  with 
its  collectors  opened  perpendicular  to  the  flow. 

Each  collector  contained  six  beryllium-copper  foils  which  acted  as  film.  The  foils  were 
stacked  on  top  of  each  other  in  a  unit  called  a  foil  cassette.  Each  foil  could  be  flipped  out  of  the  way, 
exposing  the  next  foil  of  the  cassette.  The  experiment  was  designed  to  expose  each  foil  for  a  period 
of  several  weeks.  Then  a  timer  would  trigger  the  end  of  exposure  for  the  previous  foil  and  expose  a 
new  foil.  This  enabled  a  temporal  and  spatial  analysis  of  the  distribution  of  particles  in  the 
interstellar  gas.  If  everything  had  gone  as  planned,  the  foils  would  be  retrieved  after  one  year  in 
space,  by  a  subsequent  shuttle  mission.  However,  the  Challenger  tragedy  delayed  the  shuttle  mission 
schedule.  As  a  result,  the  LDEF  spacecraft  orbited  for  almost  six  years.  Figure  1.1  shows  a  simplified 
outline  of  the  LDEF  spacecraft,  with  locations  of  the  IGE  collectors  sketched  in. 
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Fig.  1.L  A  simplified  sketch  of  the  lOE  experiment  on  the  LDEP  spacecraft  Each  square  segment 
represents  a  trs^  on  ^ch  all  experiments  were  attadied,  below  is  an  expanded  diagram  of  a  tray. 
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After  recovery  of  LDEF,  a  new  surprise  was  discovered.  The  initiate  relay  for  the  exposure 
timer  had  stuck.  Because  this  relay  failed,  the  foils  never  cycled  through  their  scheduled  exposure 
times.  Each  collector  exposed  only  one  foil  for  the  entire  flight  The  experiment  designers 
recalculated  a  computer  simulation  with  this  new  information.  Some  of  the  temporal  and  spatial 
information  was  lost  by  the  relay  failure  but  the  ability  to  identify  interstellar  gas  particles  was  not 

Discriminating  between  atmospheric  noble  gas  constituents  and  interstellar  noble  gases  is 
crucial  to  the  experiment  The  foils  were  placed  at  the  bottom  of  each  collector  housing.  The  housing 
acted  to  filter  out  ambient  helium  particles  because  they  were  serrated  and  baffied  to  prevent  easy 
access  to  the  foils.  Five  of  the  collectors  were  oriented  perpendicular  to  the  direction  of  motion. 
Nevertheless,  a  certain  amount  of  terrestrial  helium  was  expected  to  imbed  into  the  foils.  To  better 
estimate  this  unwanted  flux  of  helium,  one  collector,  number  seven,  was  oriented  towards  the  Earth 
solely  to  measure  background  particles.  The  Earth  would  prevent  any  interstellar  particles  from 
accessing  number  seven’s  foils.  Two  collectors,  numbers  two  and  three,  were  tilted  along  the  direction 
of  motion.  The  orientation  of  all  seven  collectors  will  be  discussed  in  chapter  3  below. 

The  foils  are  currently  undergoing  mass  spectrometer  analysis  at  the  University  of  Bern, 
Switzerland.  Special  techm'ques  are  required  for  measuring  the  minute  quantities  of  interstellar  gas. 
At  first  relatively  low  temperatures  (4S0  *C)  are  used  to  release  the  least  tightly  bound  particles. 
These  low  energy  particles  will  contain  the  greatest  proportion  of  background  contamination,  helium 
from  the  atmosphere.  At  higher  temperature  steps,  particles  which  penetrated  deeper  into  the  foils 
will  progressively  be  released.  This  analysis  will  provide  information  on  how  much  kinetic  energy  the 
interstellar  gas  particles  possessed. 

Early  results  from  Switzerland  indicate  helium  and  neon  have  been  found  on  the  foils.  The 
exp^ted  particle  concentration  pattern  has  been  suggested  by  the  early  analysis.  Intensity  gradients 
across  each  foil  have  also  been  found.  This  analysis  could  provide  important  data  on  the  synthesis 
of  light  nuclei  in  the  early  universe  and  the  beginning  processes  of  nucleosynthesis  and  galactic 
evolution  [Lind  et  oL,  1991]. 
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1.2  Objectives  of  the  Research 

The  goal  of  this  research  was  to  numerically  estimate  the  number  of  atmospheric  noble  gas 
particles  expected  to  accumulate  on  each  foil  Since  each  foil  had  a  concentration  gradient  due  to  the 
shadowing  effect  of  the  collector  housing,  this  gradient  was  also  calculated.  With  the  flux  of  '’He 
particles  calculated,  ^He  concentration  was  estimated  using  known  ratios. 

The  state  of  the  thermosphere  must  be  determined  to  calculate  the  flw  of  atmospheric  helium 
impinging  on  the  foils.  In  particular,  the  temperature  and  density  of  helium  at  low  Earth  orbit  (LEO) 
altitudes  must  be  given  in  order  to  apply  kinetic  theoiy.  Since  no  measured  data  specific  to  the  LDEF 
orbit  exists  for  these  parameters,  a  model  was  adopted,  the  Mass  Spectrometer  and  Incoherent  Scatter 
model  (MSIS)  created  by  A.E.  Hedin  and  his  associates  at  Goddard  Space  Flight  Center.  MSIS  is  the 
international  reference  atmosphere  for  the  thermosphere.  It  provides  temperature  and  density  values 
for  the  main  neutral  constituents  of  the  atmosphere. 

The  problem  was  broken  into  two  parts.  The  first  part  involved  providing  the  MSIS  model 
all  the  necessary  inputs  it  required.  To  receive  the  temperature  and  density  data,  the  model  requires 
information  for  specific  location,  altitude,  local  solar  time,  and  universal  time.  The  MSIS  input 
parameters  also  require  the  F10.7  cm  flux,  and  the  geomagnetic  Ap  indicator.  The  F10.7  cm  flux  is 
an  indicator  of  total  solar  energy  output  The  Ap  indicator  is  an  average  value  which  provides 
information  on  the  geomagnetic  field.  Both  of  these  indicators  track  the  changes  in  solar  output  as 
the  solar  cycle  progresses.  The  F10.7  cm  fluxes  and  Ap  values  were  accessed  firom  the  National 
Center  for  Atmospheric  Research  archives.  Solving  the  orbit  equation  at  a  specific  time  provided  the 
exact  location,  local  solar  condition,  and  altitude  of  LDEF.  The  calculated  values  of  helium  density 
and  temperature  were  stored  for  the  integration  process. 

The  integration,  the  second  part  of  the  problem,  used  the  prindples  of  kinetic  theoiy.  The 
atmosphere  was  assumed  to  be  in  thermal  equilibrium,  with  large  numbers  of  particles  interacting 
through  elastic  collisions.  This  allowed  the  Maxwell-Boltzmann  velocity  distribution  function  to 
represent  the  state  of  the  thermosphere,  at  low  Earth  orbit  altitudes.  This  distribution  function  is 
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defined  so  that  f(r,v,t)  is  the  number  of  molecules  at  position  r  +  dr,  with  velocity  between  v  and  v 
+  dv,  at  time  t  The  thermal  equilibrium  assumption  allows  the  use  of  the  Maxwell-Boltzmann 
velocity  distribution  function: 

-  (i»i/2nJk7)*^exp(-mv*/2*7)  (1.2.1) 

Where  m  is  the  mass  of  the  particle,  T  is  the  temperature,  k  is  Boltzmann’s  constant,  and: 

V  is  the  velocity  of  the  particle.  By  integrating  the  Maxwell-Boltzmann  velocity  distribution  function 
over  all  possible  velocities,  every  particle  is  represented  or  accounted  for. 

M 

At  this  point,  MSIS  gives  the  number  of  helium  particles  per  unit  volume  and  the  distribution 
function  gives  the  probability  of  finding  that  particle  with  a  specific  velocity.  This  is  all  that  is  required 
to  determine  what  firaction  of  the  helium  molecules  are  able  to  hit  the  foil.  An  experimentally 
determined  function  was  included  to  find  the  probability  that  a  particle  hitting  the  foil  would  actually 
adhere  in  the  lattice.  Various  differences  in  look  angle  were  modelled  by  changing  the  limits  of 
integration.  This  permitted  the  calculation  of  the  concentration  gradient  across  each  foil.  The  limits 
of  integration  were  modified  to  represent  the  seven  difierent  fields  of  view  and  the  difierences  in  the 
geometry  between  collectors. 

Since  this  is  a  probabalistic  analysis  special  care  must  be  taken  to  account  for  sources  of  error. 


Uncertainties  will  arise  due  to  assumptions  in  the  derivation  of  the  Maxwell-Boltzmann  equation,  the 
MSIS  model  results,  and  in  the  function  for  the  probability  of  the  particle  imbedding  in  the  foil  based 
on  it’s  momentum.  The  error  analysis  will  be  detailed  in  chapter  4. 
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CHAPTER  2 

THE  LOW  EARTH  ORBIT  ENVIRONMENT 

The  LDEF  spacecraft  orbited  the  Earth  for  almost  six  years  ifl  the  portion  of  the  atmosphere 
known  as  the  thermosphere.  This  region  responds  to  changes  in  solar  activity.  The  discussion  below 
will  examine  how  solar  activity  changes  the  density,  temperature,  and  composition  of  the 
thermosphere.  Then  it  will  examine  the  MSIS  model  and  finally,  turn  to  the  geometry  of  the 
experiment 


2.1  The  Thermosphere 

The  Interstellar  Gas  Experiment  was  launched  into  low  Earth  orbit  (LEO)  in  the  bay  of  space 
shuttle  Challenger,  in  April,  1984.  On  April  7,  the  crew  deployed  LDEF  into  a  stable  orbit  at  475 
km  altitude.  When  LDEF  was  retrieved  in  1990,  its  altitude  had  decreased  to  320  km.  NASA 
captured  it  perhaps  a  month  before  it  would  have  re-entered  in  a  free  fall  meltdown;  therefore,  this 
description  of  the  thermosphere  will  focus  on  the  LDEF  orbit  In  particular,  it  will  concentrate  on 
the  300  to  500  km  altitude  range,  at  low  latitudes,  during  the  time  April,  1984  to  January,  1990. 

Neutral  gases  outnumber  the  charged  particles  by  six  orders  of  magnitude  at  these  heights. 
This  tenuous  atmosphere  is  in  thermal  equilibrium.  Temperatures  and  densities  respond  to  the  daily 
effects  of  the  sun,  but  slowly,  with  24  hour  cycles.  Likewise,  there  are  enough  collisions  to  keep 
temperature  and  density  spatially  continuous  and  slowly  changing.  Thermospheric  winds  rarely  exceed 
one  percent  of  the  virtual  wind  produced  by  the  spacecraft  hurtling  through  a  calm  atmosphere. 

In  the  spring  of  1984,  when  LDEF  was  launched,  the  sun’s  magnetic  activity  was  declining  as 
the  solar  cycle  waned.  From  1985  through  the  end  of  1987  key  indicators  of  solar  output  remained 
consistently  low.  This  is  called  a  "quiet  sun,"  or  solar  minimum  pericxL  The  F10.7  cm  wavelength  is 


avfl  dally  FI 0.7  cm  flux  avg  81  day  FI 0.7a  cm  flux 
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accepted  by  aeronomers  as  an  indicator  of  total  solar  extreme  ultraviolet  flux  (ELFV).  MSIS  uses  a 
daily  F10.7  cm  value  and  an  81  day  average  value  as  inputs.  Figure  2.1  presents  these  two  indicators, 
averaged  over  the  entire  month,  for  the  duration  of  the  spacecraft’s  flight  The  figure  shows  the  rising 
solar  flux  beginning  in  month  45,  January,  1988. 

The  Earth’s  magnetic  field  responds  to  changes  in  solar  output  also.  The  geomagnetic  field 
creates  large  atmospheric  disturbances  during  magnetic  storms,  particularly  at  high  latitudes.  MSIS 
models  the  neutral  components.  The  charged  particles  energize  the  neutrals  during  these  storms.  The 
geomagnetic  indicator  needed  for  the  model  is  the  Ap  index,  a  global  average  of  magnetogram 
variances.  The  Ap  index  begins  as  the  Kp  indoc,  a  logarithmic  value  taken  every  three  hours  at  12 
stations  around  the  globe.  It  is  then  converted  to  the  ap  index,  which  is  a  linear  value.  The  Ap  index 
is  the  result  of  a  daily  average  of  the  three  hour  ap  values  \Knecht  and  Shuman,  1985].  figure  2.2 
graphs  the  monthly  average  of  the  Ap  values  during  LDEFs  flight  The  peak  Ap  average  value  of 
43  occurred  in  stormy  March,  1989. 

Since  the  ambient  LEO  atmosphere  is  sensitive  to  solar  i^e  dtanges,  it’s  clear  that  daily 
sunrise  and  sunset  will  produce  profound  changes  in  the  state  of  the  gas.  At  sunrise,  the  temperature 
at  altitude  has  reached  its  nadir.  As  solar  energy  heats  the  gas  through  EUV  absorption,  the 
temperature  rises  about  50  ”C  per  hour,  peaking  just  prior  to  sunset  This  heating  produces  a 
corresponding  expansion  of  the  gas  at  local  solar  afternoon.  This  diurnal  bulge  of  high  pressure 
remains  at  the  subsolar  point  Observers  on  the  Earth  measure  a  diurnal  tide  ‘pulled*  by  the  sun. 
I^gure  23  illustrates  these  diurnal  variations  for  one  full  orbit  An  interesting  effect  can  be  seen  in 
this  figure.  The  total  atmospheric  density  peaks  in  late  afternoon  as  the  temperature  peaks,  but 
helium  density  often  peaks  toward  sunrise. 

The  plotted  orbit  has  an  inclination  of  28  degrees  with  an  altitude  of  400  kilometers.  The 
relevant  solar  indices  supplied  to  MSIS  for  August  1989  were  an  Ap  of  20  and  F10.7  cm  daily  flux  of 
223.  The  81  day  average  F10.7  cm  flux  was  219.  These  high  values  represent  Che  onset  of  solar  <ycle 
22  which  would  peak  in  1991.  From  January  1989  on,  LDEF  orbited  in  a  thermosphere  agitated  by 
an  ‘active*  sun.  Typical  F10.7  cm  fluxes  were  near  200,  and  Ap  indices  were  near  20  for  the 
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month  (Apr  *84  to  Jan  *90) 


Fig.  2J  Mean  monthly  Ap  geomagnetic  index  for  the  entire  70-month  flight  of  LDEF. 
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remainder  of  the  LDEFs  journey. 

These  two  effects,  the  daily  rise  and  fall  in  temperature  associated  with  the  Earth’s  rotation, 
and  the  long  term  influence  of  steadily  inaeasing  solar  activity,  cause  the  largest  change  in  the  state 
of  the  atmosphere  at  LEO  altitudes.  Another  modulated  change  occurs  because  of  the  Earth’s 
revolution  about  the  sun.  There  is  a  bulge  of  higher  densities  at  constant  altitude  due  to  the 
equinoxes.  Figure  Z4  illustrates  this  effect  for  1989.  Figure  2.5  illustrates  the  same  effect  during  solar 
quiet  time  in  1986.  A  comparison  of  values  between  1986  and  1989  shows  that  the  active  sun 
produces  a  thermosphere  which  is  warmer  by  400  *C  and  richer  in  helium  density  by  a  factor  of  two. 

Another  important  effect  became  significant  in  1989  also.  The  denser  atmosphere  introduced 
larger  drag  and  LDEF  began  its  descent  in  earnest  Figure  2.6  graphs  the  decrease  in  altitude  for  the 
70-month  flight  Placed  in  orbit  at  475  km,  altitude  decreased  to  320  km  by  January,  1990.  The  slope 
of  the  curve  indicates  how  little  time  NASA  had  left  to  recapture  LDEF.  As  the  spacecraft  fell,  it 
also  accelerated.  Its  velocity  increased  from  7.6  km/sec  to  7.7  km/sec 

The  remaining  considerations  cause  smaller  changes  to  occur  in  the  state  of  the  thermosphere. 
The  MSIS  model  requires  specific  inputs  for  location,  year,  date,  universal  time,  and  local  solar  time. 
Each  piece  of  information  alters  the  computed  values  for  temperature  and  helium  density.  For 
example,  if  all  other  inputs  are  held  constant,  with  the  date  changing  by  five  days,  MSIS  returns  values 
for  temperature  and  density  which  differ  because  of  slight  changes  due  to  the  Earth’s  new  position 
in  its  orbit  about  the  sun.  The  MSIS  model’s  response  to  variations  in  the  input  parameters  will  be 
examined  in  chapter  4. 

Knowing  exactly  where  LDEF  was,  at  any  given  time,  could  be  calculated  using  the  Keplerian 
orbital  elements  found  in  table  2.1.  Notice  the  eccentricity  is  approximately  zero.  This  implies  a 
nearly  circular  orbit  The  orbit  was  assumed  to  be  circular  for  this  calculation. 

LDEFs  orbital  inclination  of  28  degrees  kept  it  at  low  latitudes.  The  thermosphere  behaves 
differently  in  the  tropics  compared  with  the  poles.  Strong  geomagnetic  storms  have  a  large  effect  on 
polar  wind  fields.  Ions,  driven  along  magnetic  field  lines,  produce  enough  drag  on  the  neutrals  to 
alter  circulation.  The  influx  of  storm  energy  dissipates  gradually  by  heating  the  neutrals.  This 


Avfl  He  den  (m'*)  avg  mean  Temp  (K) 
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Table  2.1  Kepleriao  Orbital  Elements  of  LDEF. 


Orbital  Element  (units) 

At  Launch 

At 

Docking 

Semimajor  Axis  (km) 

6861 

6715 

Eccentricity 

0.00141 

0.00152 

Inclination  (degrees) 

28.3709 

28.4826 

Argument  of  Perigee 

53.2328 

7.2505 

Rt.  Ascension  Ascnd.  Node 

228.409 

147.934 

True  Anomaly 

345.352 

0.2738 

Altitude  (km) 

474 

325 

Precession  Rate  (deg/day) 

-6.802 

-7.281 

produces  a  region  of  higher  pressure  which  affects  low  latitu :  ^  u  Jiuity  and  temperature  in  due  time. 

The  complexity  of  this  situation  can  itardiy  be  overstated.  Consider  the  following:  the  bulk 
of  the  corpuscular  and  joule  heating  v^ccurs  at  the  poles,  but  the  subsolar,  high-energy  heating  and 
direct  energy  injection  by  the  ring  current  mechanism  also  produce  a  low  latitude  density  disturbance. 
There  is  a  response  time  for  each  effect  Topically,  it  takes  four  to  six  hours  for  polar  density  waves 
to  reach  low  latitudes.  The  magnitudes  of  the  density  disturbances  are  also  damped  with  time  [Prolss, 
1982]. 

It  has  also  been  shown  that  different  chemical  constituents  respond  to  storms  differently. 
Helium  and  atomic  oxygen  densities  decrease  in  the  auroral  zone  and  increase  at  low  latitudes. 
Nitrogen  increases  quickly  in  direct  relation  to  geo-magnetic  activity.  Height  variations  also  arise,  and 
atomic  o^gen  is  depleted  at  low  altitudes,  but  increases  at  higher  altitudes  where  the  effect  of 
increased  temperature  is  more  important  than  transport  effects  [Engebretson  and  Afauersberger,  1983; 
Hedin  et  aL,  1977;  Hedin  ct  aL,  1981]. 
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The  climatological  circulatioo  of  the  thermosphere  beh:}ves  like  a  Hadl^  cell  during  quiet 
sun  conditions.  This  cellular  flow  features  a  low  latitude  region  of  rising  warm  air,  a  high  altitude 
poleward  flow,  sinking  at  the  poles,  and  finally  a  low  altitude  southward  flow.  During  active  periods, 
storms  produce  enough  heat  to  create  a  seomd  cell  pattern,  the  opposite  of  the  Hadley  circulation. 
In  the  tropics,  average  storm  circulations  are  damped  enough  to  show  little  effect  on  the  Hadlqr  flow. 
This  flow  is  generally  two  orders  of  magnitude  less  than  LDEFs  orbital  velodty  [Killeen  and  Roble, 
1988].  The  low  latitude  effect  called  the  equatorial  electrojet  is  capable  of  generating  zonal  winds  of 
over  150  m/s.  This  strong  current  occurs  below  120  km  altitude  [Fqer,  1991].  Its  effect  at  LDEFs 
altitude  is  not  well  understood.  Both  the  electrojet  and  the  Hadley  circulation  were  assumed  to  be 
negligible  for  this  project 

Thermospheric  research  continues  to  find  new  relationships  about  how  incoming  energy 
effects  the  density,  composition,  and  temperature  of  the  gas.  The  atmosphere  responds  both  spatially 
and  temporalfy.  This  complex  behavior  has  been  modelled  by  the  MSIS  programmers.  Aeronomers 
have  relied  on  this  semi-empirical  model  for  almost  twenty  years.  A  close  look  at  MSIS  results  will 
help  to  better  understand  LDEFs  orbital  environment  The  following  section  examines  the  MSIS 
model,  including  a  brief  look  at  how  it  works. 

2.2  The  MSIS  Thermospheric  Model 

What  is  a  "model  atmosphere"?  It  can  consist  of  tables,  graphs,  or  software  which  show 
distributions  of  pressure,  temperature,  and  gas  concentration.  Neutral  models  are  based  on 
experimental  data  and  fitted  to  the  data  using  least  square’s  best  fit  techniques.  Some  neutral  models 
in  use  today  include  the  U.S.  Standard  Atmosphere,  the  Jacchia  Reference  Atmosphere,  and  the 
Atmospheric  Handbook.  Continuous  fields  may  now  be  interpolated  from  known  data  points 
according  to  the  principles  of  an  ideal  gas  in  hydrostatic  equilibrium.  Since  the  thermosphere 
responds  to  changes  in  solar  activity  and  geomagnetic  activity,  an  effective  model  must  include  the 
option  to  input  the  proper  indicators  of  this  activity. 
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The  earliest  models  of  the  neutral  thermosphere  were  based  primarily  on  total  densities 
inferred  from  satellite  drag  analysis  done  by  L.  O.  Jacchia  in  the  early  1960s.  Incoherent  scatter  data 
showed  significant  differences  with  the  results  of  these  models  [Sojka,  1989].  Hedin’s  MSIS  model 
began  with  the  analysis  of  in  situ  measurements  from  a  neutral  mass  spectrometer  aboard  the  Ogo-6 
satellite.  Ogo-6  flew  from  June,  1969  until  July,  1971.  Its  perigee  was  at  400  km  with  an  apogee  of 
1100  km.  The  instrument  measured  molecular  nitrogen,  atomic  oxygen,  and  helium.  Hedin’s  1973 
model  provided  a  faithful  representation  of  thousands  of  measured  data  points.  As  a  bonus,  it  also 
revealed  many  global  distribution  features  which  were  inherent  in  the  data  but  were  never  as  apparent 
as  this  before  [Hedin  et  aL,  1974]. 

Hedin  gathered  more  results  from  satellites  and  incoherent  scatter  radars  and  incorporated 
them  into  the  model.  The  most  significant  of  these  were  from  the  Neutral  Atmospheric  Composition 
Experiment  and  the  Open  Source  Spectrometer  experiment  aboard  the  AE-C  satellite.  He  used  the 
incoherent  scatter  findings  to  improve  the  exospheric  temperature  fit  He  added  three  new 
constituents  into  the  model:  argon,  hydrogen,  and  molecular  oxygen  [HetUn  et  oL,  1977].  This  influx 
of  new  measurements  tightened  the  outputs  adherence  to  real-time  in  situ  data.  McPherson  and 
Rishbeth  used  Malvern  data  and  subjected  MSIS  to  a  test  in  1979.  Their  objective  was  to  examine 
the  temperature  portion  of  the  model  with  data  obtained  outside  the  range  from  which  it  was 
constructed.  They  found  only  minor  discrepancies  [McPherson  and  Rishbeth,  1979]. 

At  this  point  the  weakest  temperature  modelling  agreement  came  at  high  latitudes  for  high 
magnetic  activity.  Another  weak  area  was  in  composition  measurements  at  altitudes  below  140  km. 
The  biggest  upgrade  in  MSIS-86  came  from  including  new  data  from  the  Dynamics  Explorer  satellite, 
DE-B.  This  provided  strong  coverage  of  the  polar  regions.  Terms  were  added  to  better  represent 
seasonal  variations  near  the  poles  under  both  quiet  and  magnetically  disturbed  conditions.  N,  atomic 
nitrogen,  was  also  added  to  the  menu.  The  MSIS-86  model  is  the  accepted  CXJSPAR  International 
Reference  Atmosphere  (CIRA)  for  the  thermosphere  [Hedin,  1987, 1991]. 

The  most  recent  version,  MSlSE-90,  includes  the  mesosphere  and  lower  stratosphere.  It 
provides  aeronomers  with  "a  single  analytic  model  for  calculating  temperature  and  density  profiles 
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representative  of  the  climatological  average  for  many  geophysical  conditions*  [Hedin,  pg.  1159, 1991]. 
For  the  helium  calculations  below,  both  the  1986  and  90  editions  were  used.  I  asked  Dr.  Hedin  about 
this.  He  replied,  *In  general,  MSISB-90  is  virtually  the  same  as  MSIS-86  above  120  km*  \personat 
correspondence,  1991].  Hedin  concludes.  The  model  represents  current  knowledge  of  the 
climatological  average  reasonably  well*  {personal  correspondence,  1991].  It  certainly  has  increased  the 
level  of  complexity  since  the  original  satellite  drag  models  of  the  1960’s. 

2  J  The  Physics  and  Math  of  the  MSIS  Model 

Hedin’s  group  faced  the  task  of  reproducing  the  temperature  and  composition  profile  of  the 
neutral  thermosphere  for  any  location,  at  any  time,  under  any  solar  and  geomagnetic  condition. 
Armed  with  the  laws  of  classical  mechanics,  and  assuming  an  ideal  gas  in  hydrostatic  equilibrium,  the 
Goddard  team  proceeded  to  develop  the  model’s  code. 

Th^  knew  the  largest  variations  from  climatological  mean  values  were  caused  by  solar 
activi^.  The  sunspot  record  has  been  documented  since  the  1700s.  The  periodic  fluctuations  of  the 
solar  cycle  also  occur  in  the  extreme  ultraviolet  (EUV)  wavelengths.  EUV  radiation  absorbed  in  the 
thermosphere  creates  the  bulk  of  the  heating  which  drives  neutral  temperatures  from  200  K  to  well 
over  1000  K  every  day.  A  review  glance  at  figure  Z1  shows  how  the  F10.7  cm  flux  varied  across  the 
six  year  orbit  of  LDEF.  A  preview  glance  at  figure  3.1,  a  diagram  of  MSIS  computed  temperature  and 
density  values,  shows  the  same  trend.  Qearly,  the  solar  EUV  output  forces  the  greatest  response  by 
the  thermosphere. 

The  spherical  shape  of  the  Earth  causes  differential  heating  with  respect  to  latitude.  The 
Earth’s  rotation  and  revolution  create  large  variations  in  insolation  with  respect  to  local  time,  and 
time  of  the  year.  The  atmosphere  reacts  to  strong  heating  at  the  subsolar  point  by  expanding  and 
creating  a  diurnal  tidal  bulge.  This  creates  a  pressure  forcing  which  drives  neutral  winds.  The  model 
doesn’t  attempt  to  reproduce  the  dynamics  of  the  thermosphere  but  it  must  consider  these  winds 
because  they  cause  tidal  changes  in  the  density  at  any  location.  Another  global  feature  of  the  low 
latitude  neutral  atmosphere  is  the  midnight  pressure  bulge  formed  by  thermal  tides.  This  bulge 


19 


rotates  at  the  anti-solar  meridian  [Hemro  and  Meniwether,  1980].  Hedin  sorted  the  data  into  three 
hour  blocks  to  preserve  the  (^clic  tides  with  harmonics  of  24-,  12-,  and  8-hour  periods. 

The  complexities  of  modelling  the  climate  of  the  thermosphere  are  imposing.  The  code 
development  began  with  first  principles,  so  Hedin  and  his  associates  began  with  the  ideal  gas  law: 

p  ■  nkT 

p  is  the  pressure,  n  is  the  number  density,  k  is  Boltzmann’s  constant,  and  T  is  the  temperature. 
Another  fundamental  equation  needed  was  the  hydrostatic  equation: 

Here  h  is  the  height,  m  is  the  particle’s  mass,  and  g  is  the  gravitational  acceleration.  These  equations 
can  be  combined  to  produce: 

-i.  ^ 

p  dk  kT  ‘  H 

with  H  in  meters  defined  as  the  pressure  scale  height,  la  the  thermosphere,  for  a  given  constituent, 
both  g  and  T  will  be  functions  of  height  The  foundation  of  all  neutral  density  models  rests  on  the 
faa  that  equation  23.1  is  consistent  at  any  level  where  diffusive  separation  exists  and  where  the 
distribution  is  determined  only  by  the  verticle  pressure  gradient  balandng  against  gravity.  Below  100 
km  altitude,  equation  2.3.1  holds  if  a  weighted  average  mass  is  used  to  represent  the  masses  of  all  the 
atmospheric  constituents. 

Integrating  equation  23.1  with  height  gives  the  variation  of  pressure  with  height  above  any 
chosen  base  level,  h’,  at  which  pressure  is  p»: 

(2J.2) 

a 

When  this  expression  is  iterated  to  account  for  changes  in  gravity  and  temperature  with 
height,  it  gives  reasonable  results  for  number  densities  of  constituents.  When  number  densities  are 
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known,  along  with  a  base  level  temperature,  equation  2.3.2  can  be  solved  for  the  unknown 
temperature.  A  note  of  caution  here:  there  is  no  guaranty  that  diffusive  equilibrium  holds.  Gravity 
waves  can  induce  turbulence  which  will  often  bring  low  level  density  values  to  higher  regions.  In 
general,  the  results  of  this  model  will  be  ciiroatologically  reasonable,  not  exact 

The  task  of  keeping  the  results  close  to  real  values  involved  fitting  curves  to  thousands  of  in 
situ  data  points  gathered  from  satellites  and  measurements  ft-om  incoherent  scatter  radars.  The  model 
uses  an  exponential  temperature  profile  as  a  function  of  geopotential  height  for  the  upper 
thermosphere  and  an  inverse  polynomial  in  geopotential  height  for  the  lower  thermosphere.  These 
temperature  profiles  allow  exact  integration  of  the  hydrostatic  equation  for  a  constant  mass  to 
determine  the  diffusive  equilibrium  densi^  profile  for  each  important  species.  The  exospheric 
temperature  and  density  of  each  species  at  120  km  are  expressed  as  a  function  of  local  time,  latitude, 
longitude,  universal  time,  F10.7,  and  Ap  using  spherical  harmonics  for  latitude  versus  local  time  or 
longitude  variations  (up  to  degree  4)  with  Fourier  series  in  time  of  day  (up  to  third  harmonic)  and 
day  of  year  (up  to  2nd  harmonic),  and  polynomials  for  variations  with  F10.7  index  and  Ap  index 
[Hedin,  1991]. 

The  merit  of  an  empirical  model  based  on  a  spherical  harmonic  analysis  is  that  these 
functions  are  mutually  orthogonal.  Th^  form  a  complete  set  that  can  in  principle  represent  any 
degree  of  complexity  in  the  data  by  a  ^tematic  increase  in  the  number  of  coefficients  used.  Since 
the  spherical  harmonics  are  approximate  eigenfunctions  in  the  thermosphere,  relatively  few  terms  are 
needed.  Analysis  in  terms  of  these  components  can  aid  in  understanding  the  origin  of  various 
thermospheric  variations  [HaUn  et  aL,  1974;  Mayr  and  Volland,  1970]. 

Figure  in  displays  MSIS  data  for  temperature  during  different  solar  conditions.  Note  how 
little  temperature  varies  with  altitude  at  LEO  heights.  Above  300  kilometers  temperature  is  nearly 
constant  This  is  called  the  exospheric  temperature.  This  figure  also  shows  how  increased  solar 
output,  indicated  by  the  F10.7  cm  parameter,  changes  the  gas  temperature.  If  the  F10.7  cm  flux  is 
iocreased  from  SO  to  250,  the  tetnpmfature  at  altitude  goes  ftom  680  *K  to  1400  *K.  Lastly,  figure  2.7 
shows  the  range  in  temperature  caused  by  local  solar  time  effects.  The  temperature  difference 
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between  local  solar  noon  and  local  midnight  is  30  percent  during  active  conditions. 

The  next  few  figures  show  the  variability  of  neutral  components  under  differing  solar 
conditions.  Figure  2.8  graphs  the  density  of  atmospheric  species  as  a  function  of  altitude  during  solar 
maximum  at  local  noon.  Likewise,  figure  2.9  shows  densities  during  solar  maximum  at  midnight. 
Figure  2.10  shows  densities  during  solar  minimum  for  local  noon.  Helium  density  hovers  near  2(10)^^ 
(m'^  during  quiet  sun.  The  density  doubles  during  active  sun.  One  important  feature  to  note  is  in 
figure  2.11,  where  the  helium  density  dominates  the  thermosphere  at  LDEF  altitudes  above  450  km, 
during  solar  quiet,  at  midnight 

The  MSIS  model  does  an  effective  job  portraying  the  make-up  of  the  thermosphere.  The 
figures  show  the  fleeting  changes  the  LEO  atmosphere  undergoes  as  it  responds  to  solar  forces.  The 
next  section  turns  to  the  physical  limits  of  the  Interstellar  Gas  Experiment 

2v4  The  Geometry  of  the  Experiment 

LDEF  carried  seven  lOE  collectors.  Each  exposed  a  thin  sheet  of  foil  to  a  different  region 
of  the  sky.  This  section  wUl  illustrate  the  regions  of  the  sky  each  collector  viewed.  It  will  also 
describe  these  collectors  in  detail. 

LDEF  can  crudely  be  described  as  a  t^linder,  with  a  height  twice  its  diameter.  The  external 
surface  of  LDEF  held  numerous  experiments.  Each  experiment  was  attached  to  a  tray  roughly  one 
square  meter  in  size.  Each  tray  was  given  a  designation  such  as  H6^  or  ElZ  The  tray  designation 
indicated  its  position  on  the  spacecraft 

A  useful  coordinate  system  to  describe  the  orbit  can  be  defined  using  the  focus  of  the  orbital 
path  and  the  direction  of  motion.  Figure  Z12  shows  how  the  orbiting  spacecraft  was  oriented.  The 
Z  axis  can  be  defined  as  the  line  connecting  the  center  of  gravity  of  the  Earth  with  the  center  of 
gravity  of  LDEF.  The  positive  Z  direction  will  be  towards  deep  space.  The  positive  X  direction  will 
be  defined  as  the  velocity  vector.  This  fixes  the  Y  coordinate  as  positive  to  the  north. 

When  LDEF  was  deployed,  the  launch  introduced  no  spin  about  the  Z  axis.  This  enabled 
each  collector  to  have  a  fixed  viewing  angle  along  the  YZ  plane.  LDEF  was  gravity  gradient  stabilized 
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Fig.  2.12.  The  orientation  of  LDEF  in  orbit  North  pole  is  visible,  as  well  as  six  of  the  collectors  of 
IGE. 

so  that  the  space  end  was  always  oriented  to  outer  space.  The  orbit  was  stable.  LDEF  orbited  about 
the  eanh  much  as  the  arrow  of  a  clock  hand  revolves  around  the  center  of  the  dodc. 

Figure  2.13  depicts  the  various  fields  of  view  of  the  seven  collectors  relative  to  LDEF  [Lind 
et  oL,  1991].  The  angles  listed  are  defined  in  the  YZ  plane.  The  zenith  angle  of  zero  degrees  is 
aligned  spaceward.  An  angle  of  +90*  is  aligned  down  the  positive  Y  axis.  The  +110*  collector  is 
shielded  from  interstellar  gas  partides  by  the  Earth  because  it  looks  below  the  local  horizon.  This 
collector  will  field  test  the  computational  results  of  this  paper.  Besides  the  +24*  and  the  +70° 
collectors,  two  collectors  point  below  the  plane  of  orbit,  one  at  -24*,  one  at  -70*. 

The  remaining  two  collectors  were  oriented  in  the  XZ  plane.  Figure  2.14  diagrams  the 
position  of  these  collectors  on  the  space  end  of  LDEF.  Collector  three  tipped  24*  into  the  virtual 
wind  created  as  LDEF  orbited  through  a  relatively  stationaiy  atmosphere.  The  second  of  this  pair 
was  tilted  24*  out  of  the  virtual  wind.  Together,  this  pair  allowed  measurement  of  atmospheric  helium 
so  that  a  differendng  method  could  be  used  to  solve  for  the  atmospheric  helium.  That  value  may  then 
be  subtracted  out  of  the  total  concentration.  Figure  Z15  shows  the  location  of  tray  F6  on  the  south 
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Fig.  2.13.  The  orientation  of  the  fields  of  view  of  the  seven  IGE  collector?  relative  to  LDER  Positive 
angles  are  to  the  north  of  the  LDEF  ground  track.  The  +110®  viewing  direction  was  capable  of 
collecting  only  atmospheric  helium. 

face  and  tray  E12  on  the  north  face  of  the  spacecraft  Table  2.2  summarizes  the  position  and  viewing 
angles  of  all  seven  IGE  collectors. 

The  individual  trays  held  the  collectors  in  place.  The  top  of  each  collector  opening  was  flush 
with  the  face  of  LDEF.  For  example,  no  part  of  the  collector  housing  is  exposed  above  the  plane  of 
the  spacecraft  The  trays  on  the  spaceward  face  of  LDEF,  H6  and  H9,  were  dimensioned  exactly  alike. 
Rgure  2.16  shows  the  dimensions  for  a  "glass'  model  of  tray  H6.  The  velocity  vector  of  the  satellite 
is  to  the  right,  as  indicated  by  the  coordinates.  Collector  three  is  tilted  24°  into  the  wind,  while 
collector  two  is  tilted  out  of  the  wind.  H9  was  aligned  90°  from  tray  H6,  so  the  virtual  wind  passed 
over  the  tray  baffle.  This  tray  baffle,  indicated  in  figure  2.16,  was  a  thin  strip  of  metal  located  at  the 
top  of  the  tray.  The  tray  baffles  restricted  the  fields  of  view  of  the  collectors  on  the  space  end  of 
LDEF.  The  background  particles  were  obstructed  from  entry  by  the  tray  baffle  over  foil  five.  Another 
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Table  2.2  Collector  Viewing  Angles 


Foil 

View  Angle 

Tray 

LDEF  Side 

1 

-70°  (S) 

F6 

South  Face 

2 

0°(Aft) 

H6 

Spaceward 

3 

0°  (Fore) 

H6 

Spaceward 

4 

+24°  (N) 

H9 

Spaceward 

5 

-24°  (S) 

H9 

Spaceward 

6 

-h70°  (N) 

E12 

North  Face 

7 

+110°  (N) 

E12 

North  Face 

look  at  figure  2,14  clarifies  the  orientation.  The  tray  itself  is  boxlike,  85  cm  long,  79  cm  wide,  and 
30  cm  deep. 

The  two  remaining  trays,  E12  and  F6,  were  also  dimensioned  alike,  each  on  opposite  sides 
ofLDEF.  Figure  2.16  illustrates  the  dimensions  of  tray  F6,  as  well.  The  virtual  wind  passed  aaoss 
this  tray  from  upper  left,  as  shown  by  the  coordinates.  These  trays  are  larger  than  those  on  the 
spaceward  side,  proportioned  124  cm  long,  92  cm  wide  and  31  cm  deep. 

The  most  important  determinant  of  the  geometry  of  the  IGE  was  the  collector  housing. 
Figure  2.17  displays  the  dimensions  of  one  of  these.  All  the  collectors  were  identical  to  the  one 
shown.  The  height  firom  the  top  baffle  to  the  foil  was  38  cm.  From  a  point  at  the  center  of  the  foil, 
the  top  baffle  limits  the  field  of  view.  The  opening  at  the  top  baffle  was  a  30  by  30  cm  square.  The 
foils  were  centered  at  the  base  of  the  housing.  The  foils  were  square,  22  by  22  cm.  This  suggests  a 
second  helpful  coordinate  system. 

The  coordinate  system  most  useful  for  the  calculation  of  atmospheric  helium  accretion  was 
founded  at  the  midpoint  of  the  foil.  The  positive  X  axis  was  defined  as  the  velocity  vector  of  the 
satellite  as  before.  The  Y  axis  was  drawn  at  right  angles  to  X  and  flush  with  the  foil.  The  new  Z  axis 
was  positive  pointing  out  of  the  collector  housing.  The  origin  in  this  system,  (0,0)  was  the  foil 
midpoint.  The  center  of  the  baffle  opening  was  38  cm  above  the  foil  origin  and  this  point  was  exactly 
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Fig.  2.14.  The  location  of  trays  H6  and  H9  on  the  spaceward  side  of  LDEF. 

15  cm  to  the  midpoint  of  the  top  baCQe  along  the  X  or  Y  axes.  Therefore,  the  limiting  angle  for  the 
field  of  view  from  the  foil  origin  was  the  angle  who’s  tangent  was  38  over  15,  along  Loth  axes.  This 
angle  was  roughly  22  degrees  from  zenith. 

In  this  configuration,  for  a  particle  to  hit  the  foil  origin,  it  must  have  a  negative  velocity  in 
the  Z  direction,  Le.,  v^  less  than  zero.  TTie  particle  must  also  have  an  angle  of  incidence  between  plus 
and  minus  22  degrees  in  the  XZ  plane  and  the  YZ  plane.  In  three  dimensions,  this  defines  a  field  of 
view  which  is  shaped  like  a  four-sided  pyramid  as  shown  in  figure  Z13. 

This  chapter  has  outlined  the  physical  structure  of  the  experiment,  the  parameters  of  the 
thermosphere,  and  the  MS  IS  model.  The  following  chapter  covers  the  calculations  required  to 
determine  the  amount  of  atmospheric  helium  which  accrued  on  the  foils. 
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CHAPTERS 

CALCULATIONS 

Chapter  2  dealt  with  the  physical  make-up  of  the  Interstellar  Gas  Experiment  and  its 
orientation  in  three  dimensions.  The  MSIS  model  was  examined.  Many  figures  providing  expected 
ranges  of  thermospheric  density  and  temperature  were  presented.  This  chapter  will  turn  from  the 
general  to  the  specific.  Chapter  3  will  outline  how  the  orbital  motion  was  analyzed  to  get  the 
necessary  input  for  the  MSIS  model.  It  will  then  show  how  the  MSIS  figures  of  temperature  and 
helium  density  were  analyzed  to  solve  for  the  actual  accruement  of  helium  particles  on  the  foils. 

3.1  Generating  Input  for  the  MSIS  Model 

The  MSIS  model  was  the  source  of  the  data  needed  to  apply  kinetic  theory  and  acquire  a 
solution  to  this  problem.  The  ten  specific  parameters  needed  to  run  the  model  are  year,  day,  second, 
altitude,  geodetic  latitude,  geodetic  longitude,  local  apparent  solar  time,  a  three-month  average  of 
F10.7  flux,  daily  F10.7  flux  for  previous  day,  and  the  daily  magnetic  Ap  index.  The  first  half  of  the 
problem  was  to  provide  this  information  for  the  entire  six  year  flight  of  LDEF. 

The  duration  of  the  flight  was  split  into  70  month-long  intervals,  except  for  the  first  and  last 
intervals,  which  were  22  days  and  12  days,  respectively.  For  the  date  input  to  the  model,  the  midpoint 
of  the  interval  was  chosen.  This  generally  meant  the  15th  of  each  month.  The  seconds  parameter 
was  chosen  as  midnight.  The  time  input  for  MSIS  is  now  fixed  for  the  seventy  intervals. 

The  problem  became,  given  time,  where  was  LDEF  positioned?  A  review  glance  at  table  2.1 
shows  that  at  launch,  the  altitude  was  474  km  and  the  right  ascension  at  the  ascending  node  (RAAN) 
was  228.4  degrees,  with  a  precession  rate  of  -6.8  degrees  per  day.  The  altitude  of  the  satellite  was 
monitored  closely  during  1989.  At  the  time,  no  official  recovery  plans  had  been  set.  LDEF  average 
orbital  decay  rates  were  published  on  a  monthly  basis.  This  data  was  interpolated  to  provide  altitude 
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inputs  for  MSIS  during  each  interval.  Using  similar  information  on  the  changing  precession  rate,  a 
new  RAAN  was  calculated  for  the  midpoint  of  each  month. 

Surprisingly,  the  MSIS  model  is  sensitive  to  different  RAANs.  If  the  RAAN  input  is  altered 
to  step  around  the  globe,  all  other  inputs  being  equal,  the  temperature  and  density  varied  because  of 
the  RAAN’s  effect  on  local  solar  time.  By  dividing  RAAN  by  fifteen,  it  was  converted  to  hours  and 
added  to  the  time  of  day  to  calculate  the  local  apparent  solar  time.  RAAN  was  used  as  the  geodetic 
longitude  for  the  MSIS  input. 

The  RAAN  was  the  starting  longitude  for  one  "typical"  orbit.  The  geodetic  latitude  input  was 
provided  using  the  orbital  inclination  of  28.4".  Seven  of  the  ten  MSIS  inputs  are  now  set  The 
remaining  three,  the  Ap  magnetic  activity,  the  F10.7  cm  flux,  and  the  three  month  average  F10.7a  cm 
flux,  were  accessed  from  archives  at  the  National  Center  for  Atmospheric  Research. 

Rather  than  take  a  onetime,  one  spot,  surgical  average  for  the  entire  month,  I  calculated  a 
mean  taken  around  the  globe.  This  would  preserve  some  of  the  spatial  variance  in  evidence  in  figure 
2.3.  I  modified  an  existing  circular  orbit  code  [Thompson  and  Raitt,  1989].  The  code  updated  the 
seconds  input,  the  LST  input,  the  latitude,  and  the  longitude  eveiy  thirty  seconds  for  a  ninety  minute 
orbit.  MSIS  was  called  180  times,  or  eveiy  two  degrees  of  longitude.  The  model  returned  three  values 
of  interest  for  each  location,  the  total  atm(»pheric  density,  the  helium  density,  and  the  temperature. 
Ail  values  returned  were  for  the  altitude  provided  during  that  month. 

To  preserve  some  of  the  temporal  variance  present  in  figure  2.3,  the  temperature  and  density 
values  were  separated  into  two  bins.  One  bin  held  local  midnight  to  noon  values,  the  other  held  local 
noon  to  midnight  values.  The  ninety  values  in  each  bin  were  then  averaged  to  produce  one  typical 
AM  temperature  and  helium  density,  and  one  typical  PM  temperature  and  helium  density.  The  results 
are  plotted  in  figure  3.1  for  the  AM  data.  PM  data  is  plotted  in  figure  3.2. 

Other  averaging  schemes  could  have  been  adopted.  The  "two-a-day"  method  was  chosen 
because  the  next  step,  the  integration,  becomes  prohibitively  slow.  For  example,  the  integration  takes 
about  24  hours  of  calculation  with  two  bins  per  day.  Halving  the  interval  to  four  per  day,  doubles  the 
computation  time. 
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3.2  An  Early  Estimate 

The  results  of  the  orbital  averages  from  MSIS  provided  the  parameters  needed  for  the  next 
step,  the  integration.  The  question  of  how  many  helium  particles  actually  stuck  to  the  foils  after 
nearly  six  years  in  orbit  could  now  be  addressed.  An  extremely  crude  range  of  solutions  was  calculated 
using  fluid  theory. 

Due  to  the  geometry  of  the  collector  housing,  no  part  of  any  foil  was  exposed  to  the  virtual 
wind  created  by  the  spacecraft’s  flight.  The  inner  housing  baflies  and  serrations  prevented  any 
impingement  of  rebounding  particles.  Therefore,  the  results  of  fluid  theory  are  the  empty  set,  zero. 

For  an  extreme  upper  limit,  I  assumed  the  foils  were  fully  exposed  perpendicular  to  the  flow. 
Figure  3.3  displays  some  average  information  regarding  the  LDEF  orbit.  Figure  3.4  shows  values  of 
key  thermospheric  parameters  at  LDEF  altitude,  averaged  over  the  entire  69-month  period.  Data 
from  figures  3.3  and  3.4  was  used  for  this  Ttutterfly  net"  over  estimation.  The  flux  of  helium  particles 
was  calculated  by  multiplying  the  orbital  velocity  times  the  He  density.  The  number  of  particles 
collected  on  the  foil  was  calculated  by  multiplying  the  flux  times  the  total  time  in  orbit  The  result 
was  roughly  4(10)^  helium  particles  per  square  meter  of  foil.  This  also  assumes  every  particle  will 
adhere  to  the  foil.  Qearly,  fluid  theory  had  to  be  abandoned,  but  the  extreme  limits  it  provided 
proved  useful. 


33  Application  of  the  Kinetic  Theory  of  Gases 

The  problem  fits  the  principles  of  the  kinetic  theory  of  gases,  pioneered  by  Ludwig  Boltzmann 
and  James  Clerk  Maxwell.  The  three  basic  assumptions  required  to  apply  this  theory  are  first,  that 
the  gas  has  a  large  number  of  spherical  particles  which  make  elastic  collisions  with  themselves.  The 
thermosphere  contains  3(10)^^  helium  particles  per  cubic  meter,  a  reasonably  large  number.  The 
second  requirement  assumes  the  helium  molecules  are  separated  by  large  distances  relative  to  their 
own  diameters.  This  condition  is  met  considering  the  mean  helium  diameter  is  little  more  than  2 
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ORBIT  FACTS 


KEY  DATES 

Start:  7  Apr  1984,  10:26  MST 
Fiitish:12  Jan  1990,  8:16  MST 

MEAN  ALTITUDE 

456  km 

MEAN  ORBIT  LENGTH 

42900  km 

MEAN  LDEF  VELOCITY 

7641  m/s 

TIME  FOR  ONE  ORBIT 

94  min 

TOTAL  NUMBER  OF  ORBITS 

32250 

TOTAL  TIME  IN  ORBIT 

69  months 
181864000  secs 


Fi{.  33.  A  summaiy  of  key  facts  of  the  LDEF  orhit 

Angstroms,  while  the  spacing  between  particles  is  roughly  7(10)'^  meters.  Finally,  the  theory  assumes 
that  the  molecules  exert  no  forces  on  each  other  except  through  collisions.  Helium  is  an  inert  gas 
which  makes  this  assumption  valid. 

Given  the  above  conditions,  kinetic  theory  concludes  that  alter  some  short  time,  the  particles 
will  come  into  a  state  of  equilibrium  with  each  other.  The  equUibrium  distribution  is  the  one  that 
remains  unchanged  by  elastic  collisions.  Maxwell  also  assumed  the  velocity  components  of  each 
panicle  were  independent  The  result  of  this  derivation  was  the  velocity  distribution  ftinction, 
equation  1.2.1,  reprinted  here.  Again,  this  distribution  function  is  defined  so  that  f(r,v,t)  is  the 
number  of  molecules  at  position  r  +  dr,  with  velocity  between  v  and  v  +  dv,  at  time  t 
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THERMOSPHERIC  FACTS 


MEAN  HELIUM  DENSITY 

Z9(10)“ 

MEAN  TEMPERATURE 

911  ®K 

PRESSURE 

3(10ylO 

MEAN  FREE  PATH 

1630  km 

MEAN  SQUARE  SPEED 

1^(10)*  mV 

MEAN  SPEED 

2190  m^ 

ROOT  MEAN  SQUARE  SPEED 

2375  m/s 


Fig.  3.4.  A  summaiy  of  typical  thermospheric  parameteis  at  orbit  altitudes,  averaged  over  the 
entire  69-month  flight  of  LOEF. 

{mf2KkTf^exp(,~mv*l2kJ)  (L2-1) 

Where  m  is  the  mass  of  the  particle,  T  is  the  temperature,  k  is  Boltzmann’s  constant,  and  v  is  the 
velocity  of  the  particle.  By  integrating  the  Maxwell-Boltzmann  velocity  distribution  function  over  all 
possible  velocities,  every  helium  particle  is  accounted  for. 

The  assumption  of  thermal  equilibrium  implies  any  position  r  of  given  altitude  is  similar  to 
the  next  It  also  means  that  changes  in  time  are  slow  enough  to  be  treated  as  a  constant  of  the 
integration.  The  problem  is  reduced  from  one  with  seven  variables  to  one  with  three  velocity 
variables.  By  transposing  the  geometry  of  the  collector  housing  into  velodty  space  with  certain  fields 
of  view,  the  stage  is  set  to  solve  for  the  number  of  particles  which  hit  the  foil. 
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ZA  A  Spherical  Coordinates  Estimation 

Spherical  coordinates  seemed  like  a  useful  approach  at  first  Some  error  would  be  introduced 
by  rendering  the  actual  geometry  into  spherical  geometry.  The  rectangular  opening  of  the  housing 
was  approximated  by  a  circular  opening  of  equal  area.  The  field  of  view  went  fi'om  a  pyramid  to  a 
cone.  The  advantage  of  the  spherical  coordinates  is  the  limits  of  integration  are  very  straight  forward. 

The  coordinate  system  was  aligned  as  described  above.  Tbe  origin  was  the  center  point  of 
the  foiL  The  Z  direction  was  perpendicular  to  the  foil,  positive  being  upward.  Tbe  X  direction  was 
oriented  positive  towards  the  direction  of  motion.  The  angle  phi  was  measured  from  the  v,  axis.  The 
angle  theta  was  measured  from  the  v,  axis. 

The  derivation  of  the  Maxwell-Boltzmann  velocity  distribution  function  into  spherical 
coordinates  is  outlined  here.  The  equation  for  a  drifting  Maxwellian  was: 

V*  -  (v,-vj*  i  ♦  v/  ^  ♦  (v,  -  I 

The  and  are  constants  due  to  the  orbital  velocity  of  the  spacecraft  Th^  were  defined: 

The  angle  psi  was  defined  as  the  angle  from  zenith,  or  space  direction,  to  the  Z  axis  normal  to  the 
foil,  in  the  XZ  plane.  Psi  equaled  zero  for  all  collecting  foils  except  foils  two  and  three,  where  psi 
equaled  -24*  and  24”,  respectively. 

The  velocity  equation  for  the  drifting  Maxwellian  was  transformed  using  the  following  set  of 


substitutions. 
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=  r  sinB  co6(|> 

Vy  =  r  sin6  sinit> 

=  r  cosB 

The  uait  direction  vectors  were  substituted  as  follows: 

i  -  sinO  cos6  f  *  cos6  cos<^  6  -  sm(t»  ^ 

S  *  sin6  sin<^  f  *  cos6  siiu^  6  +  cos4^  $ 
i  -  cos6  f  -  sin0  6 

After  much  algebra,  it  can  be  shown  that  the  expression  for  velocity  squared  in  spherical 
coordinates  is  equation  3.2.1. 

V*  ■  V*  ♦  ♦  v]^-2v^(Vj,^6cos<^+v^cos0)  (3.2.1) 

The  Maxwell-Boltzmann  velocity  distribution  function  becomes: 

Av,.0.4>)  ■  CjBxp-(vf*v^*v^-2v^(yjaa9cos^+v^cosO)lvl) 

The  thermal  velocity  is  defined  as: 

vi  -  2kTlm 

As  before,  k  is  Boltzmann’s  constant,  T  is  the  temperature,  and  m  is  the  mass  of  the  helium  particle. 
The  constant  C„  is  defined  as: 

c,  -  (1/n  yi>^ 

This  distribution  function  was  integrated  as  follows: 
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///  ’'r  sln6<fv/l0d«|» 

jt 

The  limits  of  phi  were  0  to  360*.  Theta  ranged  from  0  to  22®.  The  velocity  extended  from  0  to  20 
km^ 

The  probability  of  the  particle  adhering  to  the  foil  is  a  function  of  The  spherical 
coordinates  method  was  deserted  because  of  the  dependence  with  velocity.  This  dependence  required 
the  Maxwellian  to  have  no  bulk  drift  Then  the  limits  of  integration  would  be  dependent  on  the 
orbital  velocity.  These  limits  proved  to  be  elusive. 

Despite  this,  I  ran  the  spherical  calculation.  I  assumed  the  temperature  and  helium  density 
were  the  mean  values  of  figure  3.4.  The  probability  of  sticking  to  the  foil  was  estimated  at  100 
percent  The  flux  was  calculated  by  multiplying  the  integrand  by  the  Z  component  of  vr.  The  final 
result  was  multiplied  by  total  time  in  orbit  This  rough  estimate  yielded  a  result  of  2(10)^^  particles 
per  square  meter  of  foil.  This  result  was  expected  to  be  an  overestimate. 

3.5  The  Probability  of  a  Particle  Adhering  to  the  Foil 

The  beryllium-copper  foil’s  trapping  probability  was  tested  by  the  IGE  team  in  the  early 
1980’s.  The  results  were  directly  related  to  the  particle  velocity.  Table  3.1  contains  the  results  of  this 
test  For  the  purposes  of  the  calculation,  I  needed  a  continuous  function  that  fit  this  data  for  any 
particle  velocity. 

The  probability  data  was  represented  using  a  group  of  five  functions,  one  for  five  different 
ranges  of  velocity.  The  coefficients  were  adopted  using  a  weighted  least-squares  method.  The 
equations  are  presented  here.  For  velocities  less  than  20  km/s: 

Til  ■  (10)'“  -  5(10)-*(v-1000)) 

For  velocities  between  20  and  25  km/s: 


n,  -  v*(10)-*’  /  (I  -  1  932(10)-^(v-20000)) 
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Table  3.1  IGE  Trapping  Probability,  The  Eta  Function 


VELOCITY  ETA 
(km/s) 

00 

0.0 

05 

10 

— 

15 

— 

20 

1.0(10)-5 

25 

4.6(10)‘* 

30 

0.0030 

35 

0.0083 

40 

0.017 

45 

0.03 

50 

0.05 

55 

0.079 

60 

0.12 

65 

0.157 

70 

0.195 

75 

0.23 

80 

0.26 

85 

0.295 

90 

033 

95 

037 

100 

0.40 

For  velocities  between  25  and  55  km/s: 

tj,  -  (2  V*  /  1.16(10)*®)a-w)  ♦  (.00188  V*  !  4.1(10)V 

where  w  is  defined  as: 


w  -  1  -  4(I0)*5  (y  -  25000) 

For  velocities  between  55  and  60  km/s: 

^4  -  v*(10)-”  /  4(1  -  5(10)-’(v-55000)) 

For  velocities  greater  than  60  km/s: 

n,  -  V*  /  2.44(10)*® 

The  results  of  this  fit  are  plotted  in  figure  3.5.  The  fit  is  re  'sonable.  A  note  of  caution, 
however,  the  data  was  unfortunately  sparse  below  20  lun^.  The  fit  assumed  an  eta  value  of  zero  at 


zero  velocity.  The  next  available  data  point  was  the  value  10‘*  at  a  velocity  of  20  km/s.  The  value  of 
the  distribution  function  is  a  maxinium  at  zero,  but  its  value  at  20  km/s  is  only  lO*^  for  a  thermal 
velocity  based  at  the  mean  temperature  of  911  *JL  Thus,  this  function  was  the  greatest  source  of 
uncertainty  in  the  entire  calculation. 

3.6  The  Calculation  in  Cartesian  Coordinates 

3.6.1  The  Limits  of  Integration 

The  next  step  was  to  find  the  limits  of  integration  for  the  Maxwell-Boltzmann  velocity 
distribution  function,  using  velocity  space.  The  concentric  circles  of  figure  3.6  represent  an  ideal  gas 
in  thermal  equilibrium  with  zero  bulk  flow  velocity.  The  solution  of  the  velocity  distribution  function 
is  represented  by  concentric  circles  centered  at  the  origin.  The  diagram  is  two  dimensional  in  velocity 
space.  LDEP  cuts  through  this  distribution  with  a  velocity  of  7.6  km/s  along  the  direction.  Using 
the  coordinate  system  outlined  in  section  2.4,  where  X  is  the  direction  of  orbital  motion,  and  Z  is 
perpendicular  to  the  foil  face,  then  any  foil  except  cassettes  two  and  three,  can  be  represented  by  the 
point,  (Vj,  =  7.6,  Vy  =  0,  Vj  =  0  km/s).  This  argument  also  neglects  the  yaw  of  8*.  If  a  particle  has 
a  velocity  of  (7.6,0,>1},  it  is  capable  of  hitting  the  foil.  Alternatively,  if  any  particle  has  a  positive  v, 
component,  it  can  never  hit  the  foil. 

The  limits  in  v^  depend  on  the  same  way  that  the  distance  from  the  center  of  the  foil  to 
the  housing  baffles  depends  on  Z  over  X  So  the  field  of  view  in  the  XZ  plane  is  determined  by  the 
slope  of  the  line,  theta  one,  connecting  the  midpoint  of  the  foil,  to  the  midpoint  of  the  top  aft  baffle. 
The  equation  for  this  line  is: 

n  «  «,(v*  -  vj 

The  line  connecting  the  midpoint  of  the  foil  to  the  midpoint  of  the  lop  fore  baffle  is  given  by: 


b*i — ■ — ■ — ■ — I — ■ — ■ — ■ — I — ■ — ■ — ■ — I — ^ ^ — ■ — I — 
^  0  20000  40000  60000  80000 

Velocity  (m/s) 


100000 


Ftg.  3JS.  The  trapping  probabiU^  eta,  as  a  fonction  of  velod^.  The  boxed  data  points  are 
experimental  result^  the  diamond  points  are  calculated  estimates. 
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3j6.  Lines  repieseating  the  field  of  view  for  the  mUpoint  of  the  fofi  on  ooUecton  oriented 
normal  to  the  dire^on  of  motion. 


v,  -  «i  (V*  -  vj 

Theta  one  and  theta  two  were  aUowed  to  vary  daring  the  int^ration  pioocdore.  In  effect, 
the  slopes  changed  as  the  position  on  the  fofl  changed.  This  aUowed  the  calcolation  of  the  densi^ 
gradient  across  each  foQ.  The  two  lines  intersect  the  v,  axis  at  the  orbital  veloci^  of  LDEP,  v^ 

^j*  ■  ■’W  ■  760O0«/») 

The  region  is  dosed  by  the  line  at  the  Z  velod^  cutoff  For  most  calculations  the  greatest  v,  chosen 
was  IS  km/s.  In  general,  the  region  representing  the  field  of  view  was  a  triangular  wedge  as  shown 
in  figures  3.6  to  3.8. 
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Fig.  3.7.  Lines  representing  the  field  of  view  of  all  collectors  oriented  normal  to  the  direction  of 
motion  in  Vy,  v,  space,  neglecting  the  yaw  of  8°. 

The  situation  b  analogous  in  the  YZ  plane,  with  minor  variations.  For  the  five  collectors 
oriented  normal  to  the  X  direction  the  slope  of  the  field  of  view  was  also  38  units  up,  IS  units  over. 
These  slopes  were  defined  as  phi  one  and  phi  2.  Neglecting  a  slight  yaw  of  8*,  there  was  no  bulk  drift 
in  the  Y  direction  so  the  equations  defining  the  field  of  view  in  velocity  space  were: 

V,  -  ;  V,  -  ♦, 

The  velocity  dbtribution  function  proved  to  be  tymmetric  with  respect  to  the  Vy  direction,  as  shown 
in  figure  3.7.  That  fiia  enabled  the  integration  to  extend  from  zero  to  the  line  with  slope  phi  one, 
udiile  multiptying  the  sum  by  a  foctor  of  two.  The  result  of  the  limits  of  integration  was  a  three- 
dimensional  field  of  view  shaped  like  the  four-sided  pyramid  of  figure  2.13. 

The  limits  of  integration  were  slightly  different  for  coUecton  two  and  three.  These  collectors 
were  oriented  with  an  an^e  psi  equal  to  24*  in  the  direction  of  motion.  The  lines  representing  the 
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3A  Lines  representing  the  field  of  view  for  collector  three  on  tn^  H6. 

field  of  view  intersect  at  the  point  (v,  a  a  q,  b  These  fields  of  view  are  graphed  in 
figure  3.8. 

MJt  Cakalatlng  •  Density  Gradient  Across  Each  Foil 

The  limits  of  integration  were  altered  to  represent  different  positions  across  the  foil.  The 
arguments  of  the  section  above  all  describe  the  field  of  view  at  the  midpoint  of  the  foil  The  slopes 
were  always  38  cm  up  to  15  cm  over,  in  both  the  XZ  plane  for  theta,  and  the  YZ  plane  for  phL 

To  calculate  the  number  of  particles  which  could  see  the  foil  at  any  position  on  the  foil,  I 
altered  the  slopes  to  represent  the  limiting  an^es  for  the  new  position.  I  divided  the  foil  into  25 
squares.  From  the  midpoint  of  each  of  these  squares  I  recalculated  the  slopes.  These  new  slopes 
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became  the  limits  of  the  integration.  After  a  test  integration,  the  dau  showed  there  was  no  change 
in  particle  density  across  the  foil  in  the  Y  direction. 

This  simplified  the  calculation.  The  foil  was  divided  into  five  strips  parallel  to  the  Y  axis. 
The  slope,  phi  1,  used  for  each  strip,  was  fixed  at  -38flS.  Physically  this  represented  the  X  axis,  the 
centerline  of  the  foil  Assuming  the  origin  is  at  the  foil  midpoint,  the  points  integrated  to  determine 
the  density  gradient  aaoss  the  foil  were: 

(-8.9,0),  (-4.4,0),  (0,0),  (4.4,0),  and  (8.9,0) 

The  points  above  were  the  midpoints  of  the  foil  regions  which  I  numbered  1,  2,  3,  4,  and  5, 
respectively.  Figure  3.9  illustrates  the  changing  field  of  view  across  the  foil,  on  the  v^-v^  grid. 

3.4  J  The  lategndoa  Limits  Including  the  Yaw  of  8* 

When  LDEF  was  deployed,  the  spaceoaft  was  placed  into  orbit  rather  gingerly  because  any 
spin  or  wobble  introduced  would  take  a  whUe  to  stabilize.  The  shuttle  aew  did  an  excellent  Job  of 
launching  LDEF,  however,  the  satellite  was  twisted  around  the  Z  axis  with  an  angle  of  8*  to  the  right 
for  an  observer  facing  in  the  negative  Z  direction.  This  slight  twist  is  called  a  yaw.  The  crew  opted 
to  accept  this  yaw  rather  than  risk  introducing  a  more  serious  orbital  flaw  during  a  correaion. 

The  yaw  afiects  the  limits  of  'itegration.  The  intersection  of  the  lines  representing  the  field 
of  view  were  now  changed.  The  yaw  introduced  a  factor  of  sine  and  cosine  of  eight  degrees, 
multiplied  by  the  LDEF  orbital  velocity  of  7.6  kmA. 

For  example,  consider  foil  cassette  4,  on  tray  H9,  on  the  space  side  of  LDEF.  In  velocity 
coordinates  (▼pVy.v,),  the  old  intersection  of  the  fields  of  view  was  (7.6,0,0)  km/s.  The  new 
intersection  accounting  for  the  yaw  was  (73,1.0,0.4)  kmAs.  The  yaw  skewed  the  particle  distribution 
on  the  foils  from  lines  of  constant  density  which  ran  parallel  to  the  Y  axis  to  lines  with  a  slope  of 
eight  degrees.  This  yaw  slightly  alters  the  flux  on  all  foils. 

There  remains  one  singularity.  The  tray  baCfie  described  in  section  2.4  partially  obscured  the 
view  of  collector  five,  also  on  tray  H9.  Only  the  southwest  sector  of  the  foil  was  affected  by  this  tray 
baffle.  The  calculation  for  this  section  of  the  foil  needed  a  change  of  slope  for  theta  2,  in  the  XZ 
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vx 


A  profile  of  a  foO  collector  hoosing.  The  foil  cassette  is  visible  at  the  base.  The  dashed 
lines  represent  the  field  of  view  for  foil  region  1,  the  downwind  side  of  the  foil,  and  for  foil  region 
5,  the  upwind  side. 
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plane,  to  account  for  the  limiting  angle  of  the  tray  bafOe. 

3.6.4  The  Numerical  Solution 

All  the  elements  were  now  in  place  to  perform  the  integration.  The  numerical  method  used 
was  the  rectangle  rule  in  three-dimensional  velocity  space.  Below  is  a  description  of  the  program  that 
performed  the  integration. 

First  I  selected  a  low  limit  for  the  Z  component  of  velocity.  I  chose  this  limit  as  -15  km/s 
because  the  solution  to  the  distribution  function  at  this  range  is  an  extremely  low  1(10)'^.  Any 
contribution  to  the  final  sum  from  the  velocity  space  outside  this  range  would  be  insignificant  With 
this  limiting  Z  component  the  program  used  the  equations  presented  in  section  3.6.1  to  determine 
the  upper  and  lower  limits  of  integration  in  the  X  and  Y  components  of  velocity. 

The  program  was  segmented  to  calculate  two  sums  per  month,  an  AM  sum  using  the 
temperature  and  density  values  from  figure  3.1,  and  a  PM  sum  with  the  data  from  figure  32.  To  do 
this,  the  program  read  values  for  the  month  number,  the  seconds  per  month,  LDEP  orbiul  velocity, 
as  weU  as  the  temperature  and  helium  density  data  supplied  by  MSIS.  For  each  month,  this  data 
was  used  to  calculate  a  new  thermal  velocity,  and  a  new  normalization  constant,  C^.  The  new  orbital 
velocity  value  was  used  to  redetermine  and  v^^.  Another  loop  provided  new  slopes,  theta  1  and 
theta  2,  to  represent  differing  fields  of  view  for  five  different  points  on  the  centerline  of  each  foil. 

Next,  I  selected  a  small  increment  dv.  The  smallest  velocity  interval  diosen  was  SO  mfs.  This 
defined  a  cubic  increment  in  velociQr  space.  The  program  stepped  through  the  entire  range  of  velocity 
space  within  the  limits  of  integration.  With  each  step  the  program  tested  whether  a  particle  with  the 
specific  X,  Y,  and  Z  components  of  velodty  could  hit  the  foU.  If  the  particle  could  hit  the  foil,  the 
program  called  the  velocity  distribution  function  and  the  eta  function.  These  two  probabilities  were 
calculated  with  the  specific  velocity  which  could  hit  the  foil,  and  then  multiplied  together.  Also 
multiplied  to  this,  was  the  value  dr*  limes  the  Z  component  of  velocity.  These  partial  products  were 
summed  at  each  step. 

With  the  limits  of  integration  represented  in  km^,  the  integration  equation  looked  like: 
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RhO|,^  was  the  helium  density  in  particles  per  cubic  meter.  The  eta  function  was  described  in  section 
3J.  F(v)  was  defined  as: 


F(v)-  J(yJ/(y^v^  ”  (m/2itit7)*'*e*p(-«i'V2i7) 

The  units  for  the  flux  term  are  helium  particles  per  square  meter  per  second.  Figure  3.10  plots  the 
results.  The  error  estimate  is  based  on  uncertainties  in  MSIS  output  The  basic  source  of  the 
uncertainty  was  a  plus  or  minus  SO  *C  on  temperature,  and  plus  or  minus  30  percent  on  helium 
density  [Hedin  et  aL,  1974] 

The  final  sum  was  multiplied  by  half  of  the  seconds  per  month.  The  seconds  were  halved 
because  this  was  done  for  AM  and  PM  data.  After  the  normalization  constant,  C^,  and  the  helium 
density  were  multiplied  in,  the  AM  monthly  total  was  added  to  the  PM  monthly  total  for  a  final  result 
for  that  month.  Lastly,  each  month’s  subtotal  was  summed  together  for  a  final  result  with  units, 
helium  particles  per  square  meter  of  foil  surface.  Figure  3.11  displays  a  three-dimensional  view  of  the 
gradient  of  helium  density  across  the  foils. 

When  the  yaw  is  factored  into  the  calculation,  the  results  are  skewed  eight  degrees.  Figure 
3.12  displays  a  three  dimensional  profile  of  the  helium  density  with  respect  to  position  on  the  foil. 
These  results  were  calculated  for  foil  cassette  4. 

A  review  glance  at  figure  2.16  shows  the  tray  baffle  on  tray  H6.  A  similar  tray  baffle  on  tray 
H9  effected  the  distribution  on  foil  cassette  5,  because  the  virtual  wind  flowed  across  this  baffle.  On 
the  south  side  of  the  downwind  portion  of  foil  5,  the  tray  baffle  limited  the  field  of  view.  Hgure  3.13 
shows  a  profile  of  the  helium  density  with  respect  to  foil  position  for  foil  5.  Notice  how  this  figure 
compares  with  figure  3.12.  The  southwest  sector  of  the  foil  has  decreased  helium  density  because  the 
tray  baffle  intercepted  particles  which  were  free  to  hit  the  other  foils. 
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F^.  3.10.  The  beliam  partide  density  on  the  fofl  after  LDEFs  flight  Pint  data  point  is  on  the 
downwind  side  of  the  foil,  last  data  point  is  on  the  upwind  side  of  the  foil 
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3.7  The  Error  Function  Solution 


Dr  W.  Sctunickler  proposed  the  following  method  to  solve  the  problem.  The  method  applied 
the  error  function  defined  as; 

trAx)  -  dz 

JO 

This  is  similar  in  form  to  the  foil  calculation,  which  in  one  dimension  may  be  written  as: 

>tv)  «  ^1/nv^  dv 

If  I  redefine  v/vj^  as  u,  then  dv  =  Vjj,  du.  Substituting  these  values  into  the  equation  above  yields: 

Ay)  -  du 

Ja/*0 

The  Vjj,  factor  introduced  by  the  u-substitution  cancels  the  squared  term  inside  the  radical. 

Now  by  definition  of  the  error  function,  f(v)  becomes: 

Av)  -  1/2  [  erAPIv^  -  etAafvJ] 

As  before,  this  substitution  cancels  the  factor  of  square  root  pi.  That  finished  the  development  of  the 
error  function  method  in  the  one-dimensional  case. 

The  integration  for  the  foils  was  in  three-dimensional  velocity  space.  The  transformation  of 
the  numerical  integration  equation  of  section  3.6,  into  the  error  function  case  is  sketched  below. 

The  limits  of  the  integration  had  to  be  ^qiressed  in  terms  of  v^  Therefore,  the  error  function 
substitution  was  done  only  for  the  v,,  and  Vy  components  of  the  velocity.  The  v^  component  was  left 
to  vary  as  before  from  -15  km/s  to  0  km/s  in  increments  of  dv^.  This  method  transformed  the  flux 
equation  from: 

^  ’  P**  I-is  C'  V,  dv^^^ 


where: 
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F(y)-  -  imUnkTf^expi-mr^pkl) 

into: 

^  “  P*.^.  /"u  c^'(v*/vl)  n(v*+v^’'*  «(vp  dv^ 

Cg  is  the  normalization  constant: 

C,  -  (l/« 

The  alpha  equation  contains  the  error  function: 

The  substitution  constant  was  defined  as: 

K,  -  vl%IA  -  [2it77i»]it/4 

The  complication  to  this  transformation  was  due  to  the  eta  function  being  dependent  on 
Vy,  and  Vj.  To  provide  an  approximate  dependence,  I  used  the  mean  v,  value  of  v^  to  calculate  the 
V  needed  for  the  eta  function.  The  mean  Vy  value  was  zero,  again  neglecting  the  slight  yaw. 
Mathematically,  the  correct  equation  for  v  changed  from: 

V  ■  (v*  ♦  V*  ♦ 

into: 


V  •  (vi  +  V,*)"* 

Figure  3.14  plots  a  comparison  of  the  results  of  the  error  function  method  versus  the  straight 
forward  triple  integral  solution.  The  error  function  method  overestimates  the  earlier  result  by  roughly 
an  order  of  magnitude  at  region  1  on  the  downwind  side  of  the  foil.  The  results  compare  closely  on 
the  upwind  side  of  the  foil,  at  region  5. 


Foil  Region 


Fig.  3.14.  A  ccmparison  of  results  for  the  number  of  helium  partides  per  square  meter  of  foil '  The 
solid  line  represents  the  triple  Integral  solution,  the  dashed  line  represents  the  error  function  solution. 
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The  greatest  advantage  of  this  method  is  how  rapidly  it  calculates  out  The  triple 
integration  method  of  section  3.6.3  takes  over  24  hours  of  processing  time  because  of  the  repetitive 
numerical  calculations  at  each  increment  of  dv.  The  error  function  method  effectively  reduces  the 
problem  to  a  one-dimensional  integration.  Processing  time  is  measured  in  seconds,  not  hours. 

This  chapter  has  outlined  the  development  of  several  approaches  to  solving  the  problem, 
but  how  reliable  are  the  results?  The  next  chapter  will  look  at  inherent  uncertainties  in  the 


calculations. 
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CHAPTER  4 
RESULTS 

Chapter  3  presented  the  development  of  the  problem’s  solution.  This  chapter  will  examine 
possible  sources  of  error  in  this  solution.  Different  results  would  be  expected  if  a  different  procedure 
had  been  used  to  provide  the  inputs  for  the  MSIS  model  Next,  the  integration  would  provide 
different  results  if  the  temperature  and  density  inputs  varied.  This  chapter  will  look  at  several 
questions.  How  does  the  MSIS  model  react  to  ranges  in  the  input  parameters?  How  responsive  was 
the  integration  program  to  changes  in  temperature  and  density?  How  much  uncertainty  was 
introduced  by  the  lack  of  low  velocity  data  for  the  eta  function?  Fmally,  what  can  be  done  with  these 
results,  to  enhance  the  IGE? 


4.1  Variance  of  the  Results 

4.1.1  VarThig  the  MSIS  Input 

The  procedure  I  adopted  chose  one  orbit  at  midmonth  to  represent  the  thermosphere  for  the 
entire  month.  The  universal  time,  UT,  input  was  ahvs^  midnight  I  checked  how  the  output  from 
MSIS  changed  by  changing  only  the  UT  parameter  for  Juty  15, 1989.  The  orbit  program  ran  six  times 
for  that  day,  beginning  at  midnight  and  stepping  through  the  day,  every  four  hours.  The  values  for 
temperature  and  density  were  returned  and  averaged  across  the  entire  orbit  The  mean  MSIS  results 
for  temperature  and  helium  density  are  plotted  in  figure  4.1. 

Temperature  ranged  firom  1140  to  1152,  roughly  one  percent  Helium  densities  ranged  from 
2.8(10)*^  to  3.7(10)^*,  roughly  thirty  percent  In  this  particular  month  the  values  chosen  at  midnight 
would  underestimate  the  helium  density  and  overestimate  the  temperature.  This  was  due  to  the 
function’s  dependence  on  local  solar  time  and  the  longitude  input  of  RAAN.  At  other  months  the 
values  chosen  for  midnight  were  representative  of  the  midpoints  of  these  ranges.  There  were  also 
times  when  the  UT  input  would  underestimate  the  temperature  and  overestimate  the  density.  The 
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sampling  of  140  of  these  calculations  produced  the  overall  mean  of  figure  3.4. 

How  did  the  data  change  across  one  month?  For  this  experiment  I  altered  the  date  by  five 
day  periods,  beginning  on  the  third  of  the  month.  UT  was  set  at  midnight;  only  the  date  input 
changed.  The  program  averaged  the  data  around  one  complete  orbit  for  six  months  chosen  at 
random.  This  data  proved  similar  to  the  diurnal  variation  with  the  UT  effect  Temperature  ranges 
were  less  than  one  percent  during  solar  quiet  and  helium  densities  varied  as  much  as  thirty  percent 
A  review  glance  at  figure  25  shows  a  rapid  change  in  density  from  mid  May  to  mid  June,  1986.  The 
rapid  change  was  due  to  the  equinox  bulge  when  solar  energy  was  concentrated  at  the  equator.  When 
the  results  were  computed  during  a  solar  active  month,  temperature  ranges  increased.  See  figure  2.4. 
Note  the  100  *C  drop  in  temperature  from  mid  June  to  mid  July,  1989.  This  eight  percent  drop 
resulted  from  the  erratic  F10.7  indicator,  which  went  fi^om  247  in  June  to  188  in  July,  as  shown  in 
figure  2.1,  for  months  63  to  64.  Density  varied  onty  four  percent  during  this  interval. 

The  next  input  varied  was  the  right  ascension  at  the  ascending  node,  RAAN.  The  only  input 
altered  was  this  longitude,  which  was  stepped  eight  times  in  40*  inaements.  The  orbit  code  would 
use  the  longitude  input  with  the  universal  time  to  determine  the  local  solar  time.  The  output  would 
again  be  averaged  across  the  entire  orbit  This  effect  will  depend  on  the  diurnal  (yde  of  the 
thermosphere.  Figure  4.2  presents  the  results  of  this  test,  which  was  nm  on  an  equinox  month, 
March,  1988,  and  a  solstice  month,  December,  1988.  The  March  data  was  nearly  constant  as  the 
RAAN  was  altered.  The  December  density  ranged  from  3(10)*^  to  4(10)*^,  temperatures  varied  less 
than  two  percent 


4.1.2  Varying  the  Integration  Input 

The  local  solar  time  effect  was  very  important  Figure  23  shows  that  helium  density  can 
change  by  a  factor  of  three  during  one  orbit  Temperature  can  vary  by  thirty  percent  Since  the 
Maxwell-  Boltzmann  velocity  distribution  function  is  exponential  with  a  temperature  dependence, 
small  changes  in  temperature  can  be  magnified  by  the  distribution  function. 
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Fig.  4J.  The  variance  of  MSIS  computed  temperature  and  helium  density  with  the  longitude  of 
right  ascension  at  the  ascending  node,  RAAN. 
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Fig.  4.3.  A  comparison  of  integration  results  for  December,  1988,  for  foil  cassette  3.  The  local 
solar  time  was  segmented  into  eight  three-hour  bins,  and  two  12-bour  bins,  each  with  a  mean 
temperature  and  density. 

To  explore  the  LST  effect,  for  one  month  I  averaged  temperatures  into  eight  three-hour  bins. 
This  data  was  then  integrated  to  yield  the  helium  count  on  the  foil  The  data  from  this  trial  was 
compared  with  the  usual  two  12-hour  bins.  Figure  4.3  displays  the  results.  As  expected,  the  eight 
hour  integration  results  were  slightly  higher  because  of  the  higher  temperature  intervals.  The 
downwind  side  of  the  foil,  region  1,  showed  a  five  percent  increase,  the  upwind  side  resulted  in  a  no 
significant  increase.  This  test  was  repeated  for  eight  other  months  chosen  randomly,  all  with  similar 


results. 
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4.U  Varying  the  Eta  Function 

The  largest  source  of  uncertainty  in  the  development  of  this  solution  resulted  from  the  lack 
of  low  velocity  data  for  the  eta  function.  Hgure  4.4  shows  a  representation  of  two  Maxwellian  velocity 
distribution  functions.  The  top  plot  shows  a  distribution  calculated  with  a  temperature  of  1400  *K. 
The  lower  plot  was  calculated  with  a  temperature  of  700  *K.  The  particle  density  decreases 
exponentially  away  from  the  origin.  It  drops  off  so  rapidly  that  the  probability  of  a  particle  having 
a  velocity  exceeding  20  km/s  is  essentially  zero.  Yet,  it  is  only  above  20  km/s  where  the  eta  funaion 
data  has  sufficient  coverage  to  assure  a  high  confidence  fit  An  experiment  to  determine  the 
probability  of  adherence  at  low  velocities  would  be  very  helpful. 

To  illustrate  the  range  of  uncertainty  caused  by  the  eta  function,  I  developed  a  "low  eta*  and 
a  *hJgh  eta*  function.  The  overestimated  eta  was  designed  fitting  a  line  from  the  point  defined 
the  two  conditions,  velocity  equal  to  20  km^,  eta  equal  to  10'*,  (20,10’*),  to  the  point  defined  by 
velocity  equal  to  zero  and  eta  equal  to  zero,  (0,0).  The  rest  of  the  higher  velocity  probabilities  were 
equal  to  the  eta  function  described  in  section  33.  Figure  43  plots  this  high  eta  function.  The  radical 
change  of  slope  indicates  that  this  should  be  a  gross  overestimate  of  the  trapping  probability. 
Similarly,  the  underestimated  eta  extended  a  line  from  the  point,  velocity  equal  to  25  km/s,  eta  equal 
to  4.6(10)'^  through  the  point  (20,10'*).  This  line  produced  a  zero  probability  of  embedding  into  the 
foU  at  any  velocity  below  19.9  km/s.  This  value  seemed  too  restrictive,  indeed  the  integration  proved 
to  be  zero  particles  for  all  foils.  Another  attempted  fit  chose  the  threshold  value  to  be  at  the  point, 
velocity  equal  to  ten  km/s,  eta  equal  to  zero.  This  low  eta  fit  was  equivalent  to  the  best  eta  except 
for  the  threshold  condition. 

This  low  eta  was  then  inserted  into  the  integration  program  and  executed.  The  results 
underestimated  the  original  eta  function  by  a  factor  of  ten  for  the  five  foils  aligned  perpendicular  to 
the  virtual  wind.  When  the  high  eta  function  was  applied  in  the  integration,  the  results  overestimated 
the  original  eta  function  results  by  five  orders  of  magnitude.  These  results  are  plotted  in  figure  4.6 
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FI{.  4J.  A  plot  of  an  overestimate  or  high  eta  (solid  line),  and  an  underestimate,  low  eta,  (dashed 
line)  of  the  probability  of  a  particle  adheriog  to  the  foiL 
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for  all  foils  perpendicular  to  the  velocity  vector.  Likewise  figure  4.7  plots  results  for  cassette  2,  and 
figure  4.8  displays  results  for  cassette  3. 

The  estimation  of  ^He  is  based  on  the  ratio  of  ^He  to  ^He.  Latest  research  findings  have  this 
ratio  at  1  to  1.42(10)"^  [Anders  and  Grevesse,  1989].  This  ratio  can  be  applied  to  the  findings  above 
to  create  a  ^He  gradient  as  well. 


4.2  Summary 

The  goal  of  estimating  the  number  of  atmospheric  helium  particles  which  accrued  on  each 
foil  has  been  met  Beginning  with  the  very  wide  range  (  0  to  (10)^),  the  range  was  brought  into 
tighter  focus.  The  fit  of  the  eta  function  at  low  velocities  limits  the  confidence  of  the  final  result, 
however,  the  midpoint  of  the  foils  should  collect  approximately  3(10)’  particles  per  square  meter. 
Confidence  estimates  for  the  solution  lying  within  the  dashed  line  ranges  of  figure  4.6  through  4.8 
v/ere  placed  at  2/i.  The  answer  could  be  lower  if  the  threshold  value  for  eta  rose  from  ten  km/s.  Any 
uncertainties  due  to  MSIS,  or  averaging  methods,  were  dwarfed  by  this  eta  uncertainty. 

Nevertheless,  the  terrestrial  helium  was  considered  background  noise  for  the  collection  of 
interstellar  helium.  The  collectors  were  designed  to  filter  out  the  thermospheric  helium.  Figure  4.9 
compares  a  few  inherent  differences  between  interstellar  He  and  thermospheric  He.  These  differences 
can  be  used  to  differentiate  between  the  two,  and  isolate  the  interstellar  He.  The  key  result  of  this 
work  is  that  the  generous  upper  estimate  of  the  background  helium  density  is  still  low  enough  for  the 
detection  of  interstellar  gas  particles. 

The  gradient  across  each  foil  has  also  been  determined.  The  density  gradient  increases 
exponentially  across  the  foil  towards  the  downwind  side.  The  uncertainty  of  the  results  decreases 
towards  the  upwind  side  of  the  foil  because  the  higher  velocity  particles  which  impact  region  5  depend 
on  the  better  fitted  domain  of  the  eta  function. 

The  gradients  across  each  foil  have  a  unique  signature  predicted  by  the  model.  Dr.  J.  Sojka 
noted  that  this  trend  should  enable  accurate  subtraction  of  the  background  helium  even  if  the  absolute 
quantity  is  unknown. 
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Cassette  3 


Fig.  4  A  The  results  for  cassette  3  using  the  error  function  method.  The  dashed  lines  represent 
lii^ts  of  uncertainty  due  to  the  eta  function. 
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COMPARISON  OF  INTERSTELLAR  He  AND  ATMOSPHERIC  He 

INTERSTELLAR  ATMOSPHERIC 


Velocity 

24km/s 

0  km/s 

Relative  Velocity 

24  km/s 

7.6  km/s 

Temperature 
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Density 
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2.9(10)*^  m 
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Fig.  4.9.  A  comparisoa  of  parameters  of  interstellar  helium  and  atmospheric  helium. 


4J  Suggestions  for  Future  Work 

This  calculation  really  needs  a  better  set  of  data  for  the  probability  of  a  particle  sticking  to 
the  foil  at  low  velocities.  If  the  data  from  such  an  experiment  could  hold  uncertainties  to  less  than 
a  factor  of  two,  the  results  from  the  flux  calculation  could  be  narrowed  to  plus  or  minus  one  order 
of  magnitude. 

One  approach  which  might  determine  the  trapping  probability  at  low  particle  velocities  was 
suggested  by  Dr.  J.  Sojka.  Foil  seven  on  tray  E12  faced  below  the  local  horizon.  The  only  helium 
particles  which  could  have  accrued  on  it  were  background  from  the  thermosphere  and  particles 
precipitating  along  field  lines  from  the  magnetosphere.  Since  the  model  predicts  a  specific  gradient, 
this  would  act  as  a  fingerprint.  The  atmospheric  helium  actually  found  on  this  foil  could  be  compared 
to  the  model’s  predicted  values.  Then  the  difference  between  these  values  can  be  used  to  determine 
a  more  accurate  eta  function. 

The  problem  was  well  defined.  The  method  for  the  calculation  resulted  from  the  specific 
geometiy  of  the  seven  collectors.  Each  collector  had  its  own  field  of  view  so  the  limits  of  integration 
were  changed  to  model  the  specific  view  angles.  The  model  was  flexible  enough  to  handle  these 
changes. 
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This  method  can  be  applied  to  any  passive  collection  system  as  long  as  it  behaves  as  an  ideal 
gas  in  thermal  equilibrium.  It  could  be  applied  to  find  the  flux  of  any  neutral  species  given  by  MSIS; 
for  example,  those  shown  in  figure  2.10.  To  do  this,  the  probability  of  the  particles  adhering  to  the 
foil  data  would  have  to  be  provided. 

One  aspect  which  could  be  looked  at  is,  what  if  the  particles  were  allowed  one  bounce  prior 
to  impacting  the  foil?  This  would  take  a  new  set  of  view  angles  into  consideration,  but  the  bulk  of 
the  calculation  should  remain  the  same.  In  this  case,  the  foil’s  density  distributiop  vrould  show  a 
second  peak  on  the  upwind  side  of  the  foil. 

Another  feature  of  this  analysis  is  the  difference  in  calculating  time  between  the  error 
function  method  and  the  straight  ahead  triple  integral  calculation.  Depending  on  the  accuracy 
required,  a  quick  method  for  estimating  any  flux  of  this  type  would  be  to  do  two  or  three  slow  runs 
with  the  more  accurate  triple  integration.  Then  do  a  few  fast  runs  with  the  error  function  method. 
After  comparison  of  the  two  results,  create  some  fitting  function  which  will  transform  the  error 
function  results  to  fit  the  triple  integration  results.  'The  increase  in  processing  time  could  be  very 
practical,  with  less  loss  of  accuracy  than  if  the  error  function  results  were  not  modified. 

The  IGE  mass  spectrometer  analysis  is  currently  at  a  preliminary  stage.  The  isotopes 
measured  so  far  appear  to  be  occurring  at  expected  relative  amounts.  Two  foils  have  been  analyzed 
for  an  intensity  gradient  A  gradient  has  been  found  to  exist  and  is  in  the  correct  direction  [Lind  et 
aL,  1991].  The  gradients  predicted  by  the  model  calculation  are  important  They  will  allow 
discrimination  between  atmospheric  and  interstellar  helium.  The  critical  isotope  ratios  can  then  be 
identified  in  the  interstellar  gas.  These  key  ratios  may  reveal  secrets  of  the  early  universe  and  of 
nucleosynthesis  inside  distant  stars. 
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